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Abstract 

We analyze the two-orbital Hubbard model by means of the Composite Operator Method with the aim at studying 
the phenomenon of orbital selective Mott transition (OSMT). The model contains an interorbital interaction U' , in 
addition to the usual intraorbital one U . As warming-up approximation, we use a basis of two operators only, the 
Hubbard operators. The analysis of the density of states at the chemical potential as a function of the ratio between 
the bandwidths of the two orbitals shows the clear signature of an orbital selective Mott transition as expected. 
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The two-orbital Hubbard model has recently 
come into the limelight as a toy model for the study 
of a phenomenon that seems to interest a certain 
number of materials |l|2|3j : the orbital selective 
Mott transition. In a system composed of two elec- 
tronic species (a two orbital system) is possible 
that, under the influence of strong electronic corre- 
lations, one of the two orbital becomes insulating, 
while the other stays metallic as the whole system, 
obviously. In this manuscript, we present a prelim- 
inary study of what of this physics the Composite 
Operator Method ijy is capable to grasp within a 
simple two-pole approximation. We consider the 
following two-orbital Hubbard model: 

H=-2dY,t'^^^c\{^cZ^^)-^,Y,c\{^)c,{i) 

i,a i.a 
+ C/^Z?a« + C/'5I"l(»)"2W (1) 
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where cj, ^.(i) and Ca^cr(i) are, respectively, creation 
and annihilation electron fields with spin cr(=|,J,) 
and orbital index a(= 1,2), satisfying anticommu- 
tation canonical relations, i stands for the lattice 
vector Ri and i = (i,i). na,a{i) = 4,o-(«)ca,<T(«) is 
the particle density operator of electrons of spin a 
and orbital index a. U and U' are the intraorbital 
and interorbital Coulomb interaction, respectively, 
/i is the chemical potential, d is the dimensionality of 
the system, i'"^ the hopping integral of the a-th or- 
bital and ttij is the projection operator on nearest- 
neighbor sites. The double occupancy operator per 
orbital is defined as Da{i) = na,]{i)na,i{i)- We have 
also introduced the spinorial notation 

and c'^{i) = 'Ejat.jCaU). We will fix [/' = [/ 
according to symmetry considerations and use 
t^^) as energy unit. Following the Composite 
Operator Method prescriptions 4 in the pole- 
approximation fiavor, we introduce the projector 
operators ^a(«) = [1 - "-a(«)]ca(i) and r]a{i) = 
'^a(*)ca(*) (ti = 2) and the composite field = 
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'7i(*)7 C2(*)' ^2(0) operatorial basis in or- 
der to analyze the two-orbital system with different 
bandwidths (i*^^-' < t^^^). In this approximation the 
Fourier transform of the retarded Green's function 
G{i,j) = (R [tpii)i^Hj)]) is given by 



(m) 



(k) 



(2) 



The spectral functions (T^™^(k) and the poles 
i?m(k) can be computed [4J once the Fourier 
transform of the normalization matrix /(i,j) — 
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and of the matrix m(i,j) — ({i^^/;(i, i), i)}) 
(we here report only the non-zero entries) 



TOii(k) 



U'{n2 ~ xo) 
a(k)(l-ni+p(i))] 



ii2(k) - 2di(i)[A(i) + a(k)(p(i) - 4^0] 



TO22(k) = ({/-M)42'+C/'xo 

m33(k) = -Ai/jf +C/'(ni-xo) 
- 2dt(2) [A(2) + Q,(k)(l - n2 + 



m34(k) = 2dt(2)[A(2) + a(k)(j5(2) 
m44(k) = ({/-Ai)4^)+C/'xo 
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are known. This latter depends on six parameters. 
Five of them A('') = (Ut) - (Vavt), p'"' = 

have been fixed by algebra constrains [4 

n = 4-2(Cii+C22+C33+C44) (11) 
A^"^) = C2Vl,2a-l - C2\2a ("=1,2) (12) 
C2a-l,2a=0 (a =1,2) (13) 

where n is the total filling, Cmm' — 

and C^„/ = (^V'mV'm') ^""^ on-site and the 
nearest-neighbor-site correlation functions, re- 
spectively, na^(i) is the charge (/i = 0) and spin 
(/i = 1, 2, 3) density operator. The sixth one, 
;)^Q = i {ni(i)n2(i)) , representing the interorbital 



charge correlations, has been fixed through decou- 
pling Xo — \nin2. In Fig.[Tl we report the density 




Fig. 1. The density of states at the chemical potential 
Na{u) = 0) of both orbitals as a function of the Coulomb 
potential U for different values of the ratio R = t^^^ /t(^) at 
n = 2 and T = 0. 

of states at the chemical potential N{llj = 0) of 
both orbitals as a function of the Coulomb potential 
U for different values of the ratio R — t(^)/i(2) at 
n = 2 and T = 0. We can see that the critical value 
of the Coulomb repulsion Uc2 at which a gap opens 
in the density of states of the orbital 2, which has 
fuU bandwidth (i(2) ^ i), 

remains unchanged on 
varying the ratio R. On the contrary, Uci , the value 
of the Coulomb repulsion at which a gap opens in 
the density of states of the orbital 1, is extremely 
sensible to the value of the R and seems to obey a 
linear relationship with this latter. 

In conclusion, we have shown that the Compos- 
ite Operator Method is capable to obtain an orbital 
selective Mott transition scenario in the two-orbital 
Hubbard model already within the two-pole approx- 
imation. We need now to improve the basis in order 
to get a more realistic picture and compare our re- 
sults with experiments. 
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